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Using the most recent realistic effective interactions for nuclear matter with a smooth extrapola- 
tion to high densities including causality, we constrain the equation of state and calculate maximum 
masses of rotating neutron stars. First and second order phase transitions to, e.g., quark matter at 
high densities are included. If neutron star masses of ~ 2.2>Mq from quasi-periodic oscillations in 
low mass X-ray binaries are confirmed, a soft equation of state as well as strong phase transitions 
can be excluded in neutron star cores. 



I. INTRODUCTION 

The best determined neutron star masses are found in 
binary pulsars and all lie in the range 1.35 ± O.O4M0 
(see Thorsett and Chakrabarty 1999) except for the 
nonrelativistic pulsar PSR J1012-I-5307 of mass^ M = 
(2.1 ± O.8)M0 (van Paradijs 1998). Several X-ray bi- 
nary masses have been measured of which the heavi- 
est are Vela X-1 with M = (1.9 ± Q.2)Mq (Barziv et 
al., 1999) and Cygnus X-2 with M = (1.8 ± QA)Mq 
(Orosz & Kuulkers 1999). The recent discovery of high- 
frequency brightness oscillations in low-mass X-ray bi- 
naries provides a promising new method for determin- 
ing masses and radii of neutron stars (see Miller, Lamb, 
& Psaltis 1998). The kilohertz quasi-periodic oscilla- 
tions (QPO) occur in pairs and are most likely the or- 
bital frequencies uqpo = [1 /2tt) GM / R^^^^^ of accret- 
ing matter in Keplerian orbits around neutron stars of 
mass M and its beat frequency with the neutron star 
spin, vqpo — Vs- According to Zhang, Strohmayer, & 
Swank 1997, Kaaret, Ford, & Chen (1997) the accretion 
can for a few QPO's be tracked to its innermost stable 
orbit, Rms — &GM/c^. For slowly rotating stars the re- 
sulting mass is M ~ 2.2Mo(kHz/i/Qpo). For example, 
the maximum frequency of 1060 Hz upper QPO observed 
in 4U 1820-30 gives M ~ 2.25Mq after correcting for 
the Vs — 275 Hz neutron star rotation frequency. If the 
maximum QPO frequencies of 4U 1608-52 {vqpq = 1125 
Hz) and 4U 1636-536 {vqpo — 1228 Hz) also correspond 
to innermost stable orbits the corresponding masses are 
2.\Mq and 1.9Mq. Such large masses severely restrict 
the equation of state (EoS) for dense matter as addressed 
in the following. 

Recent models for the nucleon-nucleon interaction have 
reduced the uncertainty in the nuclear EoS allowing for 
more reliable calculations of neutron star properties, see 
Akmal, Pandharipande, & Ravenhall (1998) and Engvik 
et al. (1997). Likewise, recent realistic effective interac- 



tions for nuclear matter obeying causality at high densi- 
ties, constrain the EoS severely and thus also the max- 
imum masses of neutron stars, see Akmal, Pandhari- 
pande, & Ravenhall (1998) and Kalogera & Baym (1996). 
We will here elaborate on these analyses by incorporat- 
ing causality smoothly in the EoS for nuclear matter and 
allow for first and second order phase transitions to, e.g., 
quark matter. Finally, results are compared with ob- 
served neutron star masses and concluding remarks are 
made. 



II. THE NUCLEAR EQUATION OF STATE 

For the discussion of the gross properties of neutron 
stars we will use the optimal EoS of Akmal, Pandhari- 
pande, & Ravenhall (1998) (specifically the Argonne 
F18 -f (5w-h UIX* model- hereafter APR98), which is 
based on the most recent models for the nucleon-nucleon 
interaction with the inclusion of a parametrized three- 
body force and relativistic boost corrections. The EoS 
for nuclear matter is thus known to some accuracy for 
densities up to a few times nuclear saturation density 
no = 0.16 fm~3. We parametrize the APR98 EoS by a 
simple form for the compressional and symmetry ener- 
gies that gives a good fit around nuclear saturation den- 
sities and smoothly incorporates causality at high densi- 
ties such that the sound speed approaches the speed of 
light. This requires that the compressional part of the 
energy per nucleon is quadratic in nuclear density with a 
minimum at saturation but linear at high densities 



e = E, 



comp{n) + S{n){l - 2x) 
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Here, 



■ Up + Un is the total baryon density, x = Up/n 
the proton fraction and u = n/riQ is the ratio of the 
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baryon density to nuclear saturation density. The com- 
pressional term is in Eq. (|l|) parametrized by a sim- 
ple form which reproduces the saturation density and 
the binding energy per nucleon £q — 15.8MeV at hq 
of APR98. The "softness" parameter s ~ 0.2, which 
gave the best fit to the data of APR98 (see Heiselberg & 
Hjorth-Jensen 1999) is determined by fitting the energy 
per nucleon of APR98 up to densities of n '--^ Auq. For 
the symmetry energy term we obtain Sq = 32 MeV and 
7 = 0.6 for the best fit. The proton fraction is given by 
/3-equilibrium at a given density. 

The one unknown parameter s expresses the uncer- 
tainty in the EoS at high density and we shall vary this 
parameter within the allowed limits in the following with 
and without phase transitions to calculate mass, radius 
and density relations for neutron stars. The "softness" 
parameter s is related to the incompressibility of nuclear 
matter as Ko = 18£o/(l + s) ~ 200MeV. It agrees with 
the poorly known experimental value (Blaizot, Berger, 
Decharge, & Girod 1995), Kq ~ 180 - 250MeV which 
does not restrict it very well. From [vg/cf' — dP/d{n£), 
where P is the pressure, and the EoS of Eq. (|^), the 
causality condition Vg < c requires 



(2) 



where to„ is the mass of the nucleon. With this condition 
we have a causal EoS that reproduces the data of APR98 
at densities up to 0.6 ~ 0.7 fm~^. In contrast, the EoS 
of APR98 becomes superluminal at n « 1.1 fm~'^. For 
larger s values the EoS is softer which eventually leads to 
smaller maximum masses of neutron stars. The observed 
M ~ IAMq in binary pulsars restricts s to be less than 
0.4 — 0.5 depending on rotation as shown in calculations 
of neutron stars below. 

In Fig. |l| we plot the sound speed (ws/c)^ for various 
values of s and that resulting from the microscopic cal- 
culation of APR98 for /3-stable pn-matter. The form of 
Eq. (|l]), with the inclusion of the parameter s, provides a 
smooth extrapolation from small to large densities such 
that the sound speed Vg approaches the speed of light. 
For s = 0.0 (s = 0.1) the EoS becomes superluminal at 
densities of the order of 1 (6) fm~'^. 

The sound speed of Kalogera & Baym (1996) is also 
plotted in Fig. |l|. It jumps discontinuously to the speed 
of light at a chosen density. With this prescription they 
were able to obtain an optimum upper bound for neu- 
tron star masses and obey causality. This prescription 
was also employed by APR98, see Rhoades & Ruffini 
(1974) for further details. The EoS is thus discontin- 
uously stiffened by taking = c at densities above a 
certain value ric which, however, is lower than Ug — Buq 
where their nuclear EoS becomes superluminal. This ap- 
proach stiffens the nuclear EoS for densities ric < n < Ug 
but softens it at higher densities. Their resulting max- 
imum masses lie in the range 2.2M0^M;S2.9M0. Our 
approach however, incorporates causality by reducing the 



sound speed smoothly towards the speed of light at high 
densities. Therefore our approach will not yield an ab- 
solute upper bound on the maximum mass of a neutron 
star but gives reasonable estimates based on modern EoS 
around nuclear matter densities, causality constraints at 
high densities and a smooth extrapolation between these 
two limits (see Fig. 1). 



III. PHASE TRANSITIONS 

The physical state of matter in the interiors of neu- 
tron stars at densities above a few times normal nuclear 
matter densities is essentially unknown and many first 
and second order phase transitions have been speculated 
upon. We will specifically study the hadron to quark 
matter transition at high densities, but note that other 
transitions as, e.g., kaon and/or pion condensation or the 
presence of other baryons like hypcrons also soften the 
EoS and thus further aggravate the resulting reduction 
in maximum masses. Hyperons appear at densities typ- 
ically of the order 2no and result in a considerable soft- 
ening of the EoS, see e.g., Balberg, Lichenstadt, & Cook 
(1998). Typically, most equations of state with hyperons 
yield masses around 1.4— I.6M0. Here however, in order 
to focus on the role played by phase transitions in neu- 
tron star matter, we will assume that a phase transition 
from nucleonic to quark matter takes place at a certain 
density. We will for simplicity employ the bag model in 
our actual studies of quark phases and neutron star prop- 
erties. In the bag model the quarks are assumed to be 
confined to a finite region of space, the so-called 'bag', by 
a vacuum pressure B. Adding the Fermi pressure and in- 
teractions computed to order Ug = (7^/47r, where g is the 
QCD couphng constant, the total pressure for 3 massless 
quarks of fiavor f = u,d, s, is (see Kapusta (1988), 



ag)-B + P, + P^ 



(3) 



where Pe,/^ are the electron and muon pressure, e.g., 
Pe = /ig/127r^. A Fermi gas of quarks of fiavor i has 
density Ui — kp-fir^, due to the three color states. The 
value of the bag constant B is poorly known, and we 
present results using two representative values, B — 150 
MeVfm-3 and B = 200 MeVfm-^. We take ag = 0.4. 
However, similar results can be obtained with smaller Og 
and larger B (Madsen 1998). 

The quark and nuclear matter mixed phase described 
in Glendenning (1992) has continuous pressures and den- 
sities due to the general Gibbs criteria for two-component 
systems. There are no first order phase transitions but 
at most two second order phase transitions. Namely, at a 
lower density, where quark matter first appears in nuclear 
matter, and at a very high density (if gravitationally sta- 
ble), where all nucleons are finally dissolved into quark 
matter. This mixed phase does, however, not include 
local surface and Coulomb energies of the quark and nu- 
clear matter structures. If the interface tension between 
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quark and nuclear matter is too large, the mixed phase 
is not favored energetically due to surface and Coulomb 
energies associated with forming these structures (Heisel- 
berg, Pethick, & Staubo 1993). The neutron star will 
then have a core of pure quark matter with a mantle 
of nuclear matter surrounding it and the two phases are 
coexisting by a first order phase transition or Maxwell 
construction, see Fig. 2. For a small or moderate in- 
terface tension the quarks are confined in droplet, rod- 
and plate-like structures as found in the inner crust of 
neutron stars (Lorenz, Ravenhall, & Pethick 1993). 



IV. NEUTRON STAR PROPERTIES 

In order to obtain the mass and radius of a neutron 
star, we have solved the Tolman-Oppenheimer-Volkov 
equation with and without rotational corrections follow- 
ing the approach of Hartle (1967). The equations of 
state employed are given by the pn-matter EoS with 
s — 0.13,0.2,0.3,0.4 with nucleonic degrees of freedon 
only. In addition we have selected two representative 
values for the bag-model parameter B, namely 150 and 
200 MeVfm~^ for our discussion on eventual phase tran- 
sitions. The quark phase is linked with our pn-matter 
EoS from Eq. ([^) with s — 0.2 through either a mixed 
phase construction or a Maxwell construction, see Heisel- 
berg and Hjorth- Jensen (1999) for further details. For 
B = 150 MeVfm"^, the mixed phase begins at 0.51 fm^"^ 
and the pure quark matter phase begins at 1.89 fm""^. 
Finally, for B = 200 MeVfm"'^, the mixed phase starts 
at 0.72 fm~'^ while the pure quark phase starts at 2.11 
fm~'^. In case of a Maxwell construction, in order to link 
the pn and the quark matter EoS, we obtain for B = 150 
MeVfm^^ that the pure pn phase ends at 0.92 fm~^ and 
that the pure quark phase starts at 1.215 fm~'^, while the 
corresponding numbers for B — 200 MeVfm^^ are 1.04 
and 1.57 fm~^. 

As can be seen from Fig. || none of the equations of 
state from either the pure pn phase or with a mixed 
phase or Maxwell construction with quark degrees of free- 
dom, result in stable configurations for densities above 
~ lOriO: implying thereby that none of the stars have 
cores with a pure quark phase. The EoS with pn de- 
grees of freedom have masses M^2.2Mq when rota- 
tional corrections are accounted for. With the inclu- 
sion of the mixed phase, the total mass is reduced since 
the EoS is softer. For pure quark stars there is only 
one energy scale namely B which provides a homology 
transformation (Madsen 1998) and the maximum mass is 
Mn^ax = 2.0MQ(58MeVfm-3/B)i/2 (for a, = 0). How- 
ever, for i?x^58MeVfm~'^ a nuclear matter mantle has to 
be added and for i3^58MeVfm~'^ quark matter has lower 
energy per baryon than ^^Fe and is thus the ground state 
of strongly interacting matter. Unless the latter is the 
case, we can thus exclude the existence of 2.2 — 2.3M0 
quark stars. 



In Fig. ^ we show the mass-radius relations for the 
various equations of state. The shaded area represents 
the allowed masses and radii for i'qpo = 1060 Hz of 4U 
1820-30. Generally, 



2GM <R< 



GM 



1/3 



(4) 



where the lower limit ensures that the star is not a 
black hole, and the upper limit that the accreting mat- 
ter orbits outside the star, R < Rorb- Furthermore, 
for the matter to be outside the innermost stable orbit, 
R > Rms = QGM, requires that 



M , 



1 -I- 0.75j 
12V67rGi/QPo 

2.2Mq(1 + 0.75j) 



(5) 



kHz 

^QPO 



where j = 2ttcVsI /M"^ is a dimensionless measure of the 
angular momentum of the star with moment of inertia 
/. The upper limit in Eq. (|^) is the mass when vqpo 
corresponds to the innermost stable orbit. According to 
Zhang, Smale, Strohmayer & Swank (1998) this is the 
case for 4U 1820-30 since vqpq saturates at ^ 1060 Hz 
with increasing count rate. The corresponding neutron 
star mass is M ^ 2.2 — 2.3AfQ which leads to several in- 
teresting conclusions as seen in Fig. H Firstly, the stiffest 
EoS allowed by causality (s ~ 0.13 — 0.2) is needed. 
Secondly, rotation must be included which increases the 
maximum mass and corresponding radii by 10-15% for 
Vs ^ 300 Hz. Thirdly, a phase transition to quark mat- 
ter below densities of order 5no can be excluded, cor- 
responding to a restriction on the bag constant 5^200 
McVfm-3. 

These maximum masses are smaller than those of 
APR98 and Kalogera & Baym (1996) who, as discussed 
above, obtain upper bounds on the mass of neutron stars 
by discontinuously setting the sound speed to equal the 
speed of light above a certain density, n,.. By varying 
the density ric = 2 — > 5no the maximum mass drops 
from 2.9 2.2Mq. In our case, incorporating causal- 
ity smoothly by introducting the parameter s in Eq. (|^), 
the EoS is softened at higher densities in order to obey 
causality, and yields a maximum mass which instead is 
slightly lower than the 2.2M0 derived in APR98 for non- 
rotating stars. 

If the QPOs are not from the innermost stable or- 
bits and one finds that even accreting neutron stars have 
small masses, say like the binary pulsars, M^IAMq, this 
may indicate that heavier neutron stars are not stable. 
Therefore, the EoS is soft at high densities s^O.4 or that 
a phase transition occurs at few times nuclear matter 
densities. For the nuclear to quark matter transition this 
would require < 80 MeVfm"'' for s — 0.2. For such 
small bag parameters there is an appreciable quark and 
nuclear matter mixed phase in the neutron star interior 
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but even in these extreme cases a pure quark matter core 
is not obtained for stable neutron star configurations. 

A third QPO frequency referred to as Horisontal 
Branch Oscillations around vhbo — 20 — 50 Hz has been 
suggested to be caused by Lcnse-Thirring precession at 
the inner border of the accretion disk (Stella & Vietri 
1998) 



13.2 / i^s / t^QPo V 
50km2 M 300Hz V IkHz / 



(6) 



However, even for the stiffest EoS s ~ 0.13 — 0.2 we 
calculate moment of inertia and Lense-Thirring frequen- 
cies from Eq. (^), which are a factor ^ 4 below the 
observed vhbo thus confirming analyses of Schaab & 
Weigel (1999), Psaltis et al. (1999) and Kalogera & 
Psahis (1999). 



V. SUMMARY 

Modern nucleon-nucleon potentials have reduced the 
uncertainties in the calculated EoS. Using the most re- 
cent realistic effective interactions for nuclear matter of 
APR98 with a smooth extrapolation to high densities 
including causality, the equation of state could be con- 
strained by a "softness" parameter s which parametrizes 
the unknown stiffness of the EoS at high densities. Max- 
imum masses were calculated for rotating neutron stars 
with and without first and second order phase transitions 
to, e.g., quark matter at high densities. 

The calculated bounds for maximum masses leaves two 
natural options when compared to the observed neutron 
star masses: 

• Case I: The large masses of the neutron stars in 
QPO 4U 1820-30 (M = 2.3Mo;, PSR J 1012+5307 
(M = 2.1±0.4Mq;, Vela X-1 (M = l.giO.lM©;, 
and Cygnus X-2 (M = 1.8 ± QAMq), are con- 
firmed and complemented by other neutron stars 
with masses around ^ "ZMq . 

As a consequence, the EoS of dense nuclear mat- 
ter is severely restricted and only the stiffest EoS 
consistent with causality are allowed, i.e., softness 
parameter 0.13 < s^O.2. Furthermore, any sig- 
nificant phase transition at densities below < 5no 
can be excluded. That the radio binary pulsars all 
have masses around 1.4Mo is then probably due to 
the formation mechanism in supernovae where the 
Chandrasekhar mass for iron cores are ~ 1.5Mq. 
Neutron stars in binaries can subsequently acquire 
larger masses by accretion as X-ray binaries. 

• Case II: The heavy neutron stars prove erroneous 
by more detailed observations and only masses like 
those of binary pulsars are found. If accretion does 
not produce neutron stars heavier than ;^1.4M0, 
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this indicates that heavier neutron stars simply are 
not stable which in turn implies a soft EoS, either 
s > 0.4 or a significant phase transition must occur 
already at a few times nuclear saturation densities. 
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FIG. 1. {vs/c)'^ for /3-stable pn-matter for s = 0, 0.1, 0.2, 0.3, the results of APR98, and for the patched EoS of Kalogera & 
Baym (1996) which shows a discontinuous {va/cf ■ 
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FIG. 2. Total mass M as function of central density ric for various values of s (upper panel) and the bag parameter B (lower 
panel) for both a mixed phase and a Maxwell contracted EoS with s = 0.2 in Eq. Fjl). In addition we include also the rotational 
corrections for the pure pn-case with a = 0.2 and the mixed phase contraction for B = 200 MeVfm^'^. For the Maxwell 
construction which exhibits a first order phase transition, in the density regions where the two phases coexist, the pressure is 
constant, a fact reflected in the constant value of the neutron star mass. All results are for /9-stable matter. Note also that for 
the upper panel, the EoS for s = 0.3 and s — 0.4 start to differ from those with s — 0.13, 0.2 at densities below 0.2 fm~"^. 
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FIG. 3. Neutron star masses vs. radius for the EoS of Eq. (1) with softness s=0. 13,0.2,0.3,04, with increasing values of s from 
top to bottom for the full curves. Phase transitions decrease the maximum mass whereas rotation increases it. The shaded 
area represents the neutron star radii and masses aUowed (see text) for orbital QPO frequencies 1060 Hz of 4U 1820-30. 
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